Bulk evidence for single-gap s-wave superconductivity in the intercalated graphite 

superconductor CgYb 
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We report measurements of the in-plane electrical resistivity p and the thermal conductivity «: of 
the intercalated graphite superconductor CeYb to temperatures as low as Tc/100. When a field is 
applied along the c-axis, the residual electronic linear term no/T evolves in an exponential manner 
for Hci < H < Hc2- This activated behaviour establishes the order parameter as unambiguously 
s-wave, and rules out the possibility of multi-gap or unconventional superconductivity in this system. 
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PACS numbers: 74.70.Wz, 74.25.Fy, 74.25. Op 

Carbon is a remarkably versatile element - in its pure 
form it may exist as an electronic insulator, semiconduc- 
tor or semimetal depending on its bonding arrangement. 
When dopant atoms are introduced, superconductivrh^ 
may be added to this list, observed in graphite P, U, 
fuUerenes j2j and even diamond ^^l ■ Superconductivity in 
doped carbon was first discovered in the graphite interca- 
late compounds (GICs), materials composed of sheets of 
carbon separated by layers of intercalant atoms. The first 
of these compounds were intercalated with alkali atoms, 
and had modest transition temperatures of 0.13-0.5 K 
0. The recent discovery of T^s two orders of magni- 
tude higher than this in CeYb 5] and CgCa 'B', has 
however refocused attention on this intriguing family of 
compounds. 

The effects of the intercalant atoms in the GICs are 
two-fold: they dramatically change the electronic prop- 
erties of the host graphite lattice by both increasing the 
separation of the carbon sheets, as well as contributing 
charge carriers. This causes the two-dimensional graphite 
IT* bands to dip below the Fermi level. The graphite 
interlayer band, previously unoccupied, also crosses the 
new Fermi level, contributing three dimensional, free- 
electron like states located between the carbon sheets. 
This new interlayer band hybridises strongly with the n* 
bands, and its occupation appears to be linked with the 
occurrence of superconductivity in the GICs Q ■ 

There are still several fundamental questions remain- 
ing about superconductivity in the GICs, especially in 
CeYb and CeCa, where little experimental data exists. 
The pairing mechanism is unresolved, with speculation 
ranging from a conventional route involving the inter- 
calant phonons 0, H, 0| to superconductivity via acous- 
tic plasmons 

Early theoretical studies motivated by the alkali-metal 
GICs 1^ ^2 emphasized a two-gap model for the su- 
perconducting state. In this picture, superconductiv- 
ity arises from coupling between intercalant s electrons 
and the graphite n electrons bands, with gaps of differ- 
ent magnitudes existing on different sheets of the Fermi 



surface. Such a scenario is plausible, as there are no- 
table similarities between the GICs and MgB2 0,^3' ^ 
known multi-gap superconductor. Indeed, some aspects 
of graphite intercalate superconductivity can be under- 
stood by this two-gap phenomenology, however there is 
little direct evidence to support this picture. 

The experimental starting point for addressing these 
issues is to establish the superconducting order param- 
eter. In CeYb and CgCa, this task is complicated by 
the difficulties in materials preparation arising from the 
standard vapour transport process used for intercalation. 
This process typically yields samples which have a shell 
of fully intercalated material surrounding a core of un- 
intercalated graphite. In addition, both compounds are 
extremely sensitive to air, and their surfaces rapidly de- 
teriorate if left exposed. A recent study of penetration 
depth on CeCa suggested that the superconductiv- 
ity was s-wave, however this technique is extremely sur- 
face sensitive, and the results were dependent on surface 
treatment. With these considerations in mind, we turn 
to measurements of bulk thermodynamic properties to 
probe the superconducting state. 

The technique of thermal conductivity is ideally suited 
to the study of these materials. It is sensitive only to 
delocalized states, and in highly conductive systems such 
as CgYb the majority of the heat transport at low tem- 
peratures is provided by electrons, allowing us to easily 
separate out electronic and phononic contributions to the 
heat current. Most importantly, thermal conductivity is 
a bulk probe, only marginally affected by small concen- 
trations of impurity phases. 

In this Letter we report measurements of low temper- 
ature thermal conductivity (k) in CeYb, which we use to 
establish the nature of the superconducting order param- 
eter. The behaviour of k as the superconducting state 
is suppressed with a magnetic field shows an activated 
dependence, clear evidence of s-wave superconductivity 
with a single gap energy scale. 

Thermal transport was measured down to 60 mK in a 
dilution refrigerator using a one heater, two thermometer 
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FIG. 2: Low temperature k of CeYb as a function of applied 
field, with J|[ ab and H \\ c. The black line is the normal 
state Wiedemann-Franz law expectation estimated from the 
resistivity in H = 1 T. 



FIG. 1: In plane resistivity for CgYb with H \\ c. For H 
= T the superconducting transition is sharp with Tc — 
5.4 ± 0.4 K but broadens with applied field. For i/ = 1 T 
superconductivity is entirely suppressed, revealing a metallic 
normal state with a residual resistivity po — 5.8 fiQcm. 

steady state technique. Magnetic fields from to 1 T 
were applied parallel to the c-axis and perpendicular to 
the in-plane heat current. For measurements of k(T) at 
constant field, the sample was cooled in field from T > 
Tc to maximize homogeneity of the vortex lattice. 

Our samples were prepared by intercalating very pure 
highly oriented pyrolytic graphite (HOPG) using the 
vapour transport process described elsewhere Q. In or- 
der to obtain samples consisting mainly of intercalated 
material, we first cleaved the graphite along the a6-plane, 
and then took thin bars from the sides. The resulting 
samples were rectangular platelets of dimensions approx- 
imately 1 mm X 0.5 mm x 100 /im. Good quality electri- 
cal contacts were made using Dupont silver paint applied 
directly on the surface after cleaving, with all handling 
and mounting done in a glove box under flowing He in 
order to preserve the quality of the samples. 

Figure 1 shows the in-plane resistivity of CgYb in both 
zero field (inset) and as a function of field applied along 
the c-axis (main panel). The residual resistivity po, is 
observed to be 4.5 /illcm by extrapolating the zero field 
curve, and slightly larger than this using the 1 T curve. 
The magnetoresistance is comparatively weak in CgYb 
[l5|. only 30 % by 1 T at 4 K, compared to a factor 
of 100 increase by only 0.2 T in pure graphite at T = 
5 K . From the magnitude of po we may estimate the 
electronic mean free path £e, assuming kp ~ 0.5 A"-*^ 
we get £e — 1000 A at low temperatures in zero field. 

Figure 2 shows the dependence of the thermal conduc- 
tivity on field for H \\ c-axis. By plotting k/T vs. T 



we see that k/T is almost constant, reflecting a dom- 
inant electronic contribution. Phonons are responsible 
for the small slope, which has negligible impact of the 
field dependence of the residual linear term, the subject 
of this paper. The rise in k with H is gradual at first, 
then rapidly approaches the normal state value as Hc2 
is reached. By 1 T the sample is fully in the normal 
state, and a comparison to the electrical resistivity via 
the Wiedemann-Franz law is shown by the dashed line. 
In the normal state we see that this law is obeyed to 
within 5 %, as expected for a metal in the elastic scat- 
tering regime. 

Extrapolating the zero field data to T = in 
Fig. 2 yields a small residual linear term kq/T ~ 
0.3 mWK~^cm~^. It is tempting to attribute this as 
arising from nodal quasiparticles in the superconduct- 
ing state 01, as observed for example in d-wave super- 
conductors such as the high-Tc cuprates, where kq/T 
= 1.41 mWK-^cm-i for overdoped Tl-2201 with T^ = 
15 K . An estimate of the magnitude of the linear 
term however makes this scenario unlikely. In a nodal 
superconductor the size of kq/T is determined by ratio 
of the quasiparticle velocities parallel (^2) and perpen- 
dicular (vp) to the Fermi surface near the nodes [l3|. 
In a two dimensional d-wave superconductor with a gap 
maximum Aq and a density of n planes per unit cell of 
height c we may write 

assuming V2 0|. Usin g an average vp — 3.4 x 

10^ cm/s O, ~ 0.5 A-i]lll and Aq = 2.14kBTc 
= 1.2 meV we expect kq/T = 6.1 mWK~^cm~^ for this 
system. This is over an order of magnitude larger than 
what we measure, and on this basis we rule out an uncon- 
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FIG. 3; Thermal conductivity of CeYb at T = 300mK as a 
function of applied field. Hci and Hc2 are determined as the 
points where the conductivity deviates from a nearly flat H 
dependence at low and high fields respectively. The dotted 
line is a fit to the expected behaviour for an s-wave supercon- 
ductor between Hci < H < Hc2/2 |2l|. 



ventional order parameter in CgYb. We instead interpret 
the finite kq /T as arising from inclusions of pure graphite 
where full intercalation was not successful. The extrinsic 
origin of the linear term is supported by the fact that in 
a second sample with a similar residual resistivity (po = 
4.4 fiilcm), kq/T was seen to be half as large. 

Since thermal conductivity probes the bulk of a sam- 
ple, we may use k as a robust check of the volume fraction 
of superconducting material. Since an s-wave supercon- 
ductor has a fully gapped Fermi surface in the supercon- 
ducting state, we would expect no contribution to kq/T 
from the intercalated material. The non-intercalated 
graphite has a residual resistivity of ~ 6 fiilcm l20l | , and 
so conducts heat and charge with the same ability as bulk 
CeYb in the normal state at low temperatures. Thus, the 
ratio of kq /T in the superconducting state to that in the 
field induced normal state should, to first approxima- 
tion, yield the volume fraction of non-superconducting 
graphite inclusions. 

Performing this simple analysis yields 
{ko^h=o/T)/{ko^h=it/T) ~ 7 %, a relatively small 
fraction of the sample. These arguments demonstrate 
that good quality, essentially bulk samples of CgYb 
can be obtained by carefully selecting material from 
the edges of an intercalated platelet prepared using the 
vapour transport process. 

We now turn to analysing the field dependence of n. 
Fig. 3 shows k/T{H) with T held constant at 300 mK. 
Given the small phonon contribution evident from the 
slope of the data inFig. 2, the conductivity in the present 
figure is largely due to electrons. Starting from H = 
we see that k/T is close to zero, and essentially flat up to 
H ~ 0.025 T with a small conductivity arising from the 



non-superconducting graphite regions. For higher fields a 
sudden increase in k/T is observed, which we interpret as 
the onset of the vortex regime at H > Hd- This agrees 
reasonably well with estimates of Hd — 0.04 T at low 
temperatures from magnetization measurements 0. 

As the field is further increased, the conductivity 
evolves in an exponential manner, precisely what is ex- 
pected for transport in the mixed state of an s-wave 
superconductor. As vortices first enter the sample at 
H > Hd the only quasiparticle states at T <C Tc are 
those associated within the vortex cores When the 
vortices are far apart, these states remain localized, and 
are unable to contribute to heat transport. Increasing 
the field decreases the intervortex spacing d ^ \J<^q/B 
and the states begin to overlap, forming dispersive bands 
which yield a conductivity that grows exponentially with 
the ratio of the vortex spacing to the coherence length 
d/£_, K, oc Vi?exp(-ck-y/ Hc2/H) where a is a constant. 
This dependence is readily observed in simple s-wave su- 
perconductors such as Nb |23| for i/ci < H < i/c2/2, 
and is much different than that in nodal superconduc- 
tors, where the conductivity is observed to increase as 
^y H/ Hc2 [13 • A fit of our data to the simple s-wave 
form is shown Fig. 3, and the good agreement forms the 
central result of our work: the evolution of the electronic 
conductivity of CgYb is approximately exponential with 
applied field, providing the first verification of s-wave 
superconductivity in intercalated graphite using a bulk 
thermodynamic technique. 

At still higher fields, the conductivity rolls off and 
eventually saturates as the sample enters the normal 
state. From the data in Fig. 3 we estimate this to occur 
at ~ 0.12 T, in good agreement with estimates of Hc2 
= 0.11 T from magnetization measurements With 
Hc2{\\c) = 0.12 T we estimate £,ab ^ 525 A~ 4, which 
places CgYb in the dirty limit. This observation is con- 
sistent with the fact that the rise in conductivity with 
field is not as dramatic as in clean Nb [23, but closely 
resembles that observed in dirty limit Nb |2J| and metal 
alloy superconductors [2^ . 

In addition to confirming s-wave superconductivity, 
the activated behaviour observed in Fig. 3 rules out the 
scenario of multi-gap superconductivity originally pro- 
posed for the GIC's j^. In Fig. 4 we compare CeYb 
to other type-II superconductors by plotting the normal- 
ized value of kq/T versus applied field. The similarity 
between CeYb and the alloyed superconductor InBi with 
Tc = 4.0 K and Hc2 = 0.07 T f2j| is striking. Both 
curves are exponential with field at low H, crossing over 
to a roughly linear behaviour closer to Hc2 as exp ected 
for s-wave superconductors in the dirt y li mit |2q . The 
clean limit case observed in pure Nb [23 is shown for 
contrast. These three curves are very different from the 
behaviour of the archetypal multi-band superconductor 
MgB o 1271 ■ or the archetypal d-wave superconductor Tl- 

2201 i3r 

In the multiband scenario, gaps of different magnitudes 
are associated with the tt and a bands. Such a situation 
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FIG. 4: The normalized residual linear term kq/T of sample 
B plotted as a function of H/Hc2, with the small contribu- 
tion from graphite impurities subtracted off. For comparison, 
similar low-temperature data for the clean s-wave supercon- 
ductor Nb |23l . the dirty s-wave superconducting alloy InBi 
[25|. the multi-gap superconductor MgB2 0] and an over- 
doped sample of the d-wave superconductor Tl-2201 [3 are 
plotted alongside. 

is expected for instance when one band has strong pair- 
ing and induces superconductivity in the other by Cooper 



pair tunnelling |29j , or when electron-phonon coupUng is 
significantly different for different bands. Applying a field 
rapidly delocalizes quasiparticles states confined within 
the vortices associated with the smaller gap band, while 
those states associated with the larger gap band delocal- 
ize more slowly. This gives rise to the rapid increase in 
conductivity at low fields and relatively flat dependence 
at higher flelds ,3Q|| evident for in the MgB2 27] data 
shown in Fig. 4. Our own data thus allows us to rule 
out any sizable difference between the size of the gaps as- 
sociated with each band in CgYb, and suggests a single 
gap energy scale for the electrons as in the conventional 
s-wave scenario. 

In summary we have used bulk measurements of the 
thermal conductivity k to definitively establish s-wave 
superconductivity in CgYb, and rule out an order pa- 
rameter with nodes. The activated behaviour of kq/T 
also strongly suggests that the pairing state is isotropic, 
very similar to elementary type-II superconductors in the 
dirty limit. It will interesting to confirm these results on 
other members of the intercalate family with complemen- 
tary techniques, although it seems likely that other GICs 
will share similar superconducting properties. 
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